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t—i , ABSTRACT 

. 

| Context. IGR J16318-4848 is the prototype and one of the more extreme examples of the new class of highly obscured Galactic X-ray sources 
(**] . discovered by INTEGRAL. A monitoring campaign on this source has been carried out by XMM-Newton and INTEGRAL, consisting in three 
O l' simultaneous observations performed in February, March and August 2004. 

Aims. The long-term variability of the Compton-thick absorption and emission line complexes will be used to probe the properties of the 
circumstellar matter. 

Methods. A detailed timing and spectral analysis of the three observations is performed, along with the reanalysis of the XMM-Newton 
^ observation performed in February 2003. The results are compared with predictions from numerical radiative transfer simulations to derive the 
r* ■ parameters of the circumstellar matter. 
, ' Results. Despite the large flux dynamic range observed (almost a factor 3 between observations performed a few months apart), the source 
, remained bright (suggesting it is a persistent source) and Compton-thick (Nh > 1.2 x 10 24 citT 2 ). Large Equivalent Width (EW) emission lines 
}_i from Fe K„, Fe and Ni K a were present in all spectra. The addition of a Fe K„ Compton Shoulder improves the fits, especially in the 2004 
observations. Sporadic occurrences of rapid X-ray flux risings were observed in three of the four observations. The Fe K„ light curve followed 
the continuum almost instantaneously, suggesting that the emission lines are produced by illumination of small-scale optically-thick matter 
around the high-energy continuum source. Using the iron line EW and Compton Shoulder as diagnostic of the geometry of the matter, we 
suggest that the obscuring matter is in a flattened configuration seen almost edge-on. 
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1. Introduction 

IGR J 163 18-4848 was discovered by INTEGRAL on January 
29 th 2003, by the IBIS/ISGRI instrument, during a routine 
Galactic Plane Scan (Courvoisier et al. 2003). Re-analysis of 
archival ASCA data revealed at that position a highly absorbed 
(Nh ~ 10 24 cirT 2 ) source (Murakami et al. 2003, Revnivtsev et 
al. 2003). 

The discovery of this source triggered a XMM-Newton 
Target of Opportunity observation on February 10, 2003, and 
the subsequent analysis of the XMM-Newton data, in the 5-15 
keV energy band, confirmed the presence of a heavily absorbed 
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source, Nh ~ 2xl0 24 cm 2 , with strong emission lines (Schartel 
et al 2003). 

This complex emission could be resolved in three compo- 
nents, with centroid energies of 6.410 + 0.003 keV, 7.09 + 0.02 
keV and 7.47 + 0.02 keV, that correspond to Fe K ff , Fe Kg, and 
Ni K„ fluorescent emission lines respectively (de Plaa et al. 
2003). These preliminary results also showed that, during the 
XMM-Newton observation, IGR J 163 18-4848 was obscured 
by a Compton-thick absorber, N H = (1 -66 + 0.16) x 10 24 cirT 2 . 
A detailed analysis of this observation can be found in Matt & 
Guainazzi (2003) and Walter et al. (2003). 

Further studies at different wavelengths have also been car- 
ried out on the source. Filliatre & Chaty (2004) found the op- 
tical counterpart of IGR J16318-4848. Their studies suggest 
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Table 1. XMM-Newton monitoring campaign observations on 
IGR J16318-4848. The table shows EPIC-pn exposures times 
and source extraction region around the object for each obser- 
vation. 



Start Time - 


End Time (UT) 


T 

* exp 


Region 


2003-02-10 16:58 


2003-02-10 23:55 


23000 s 


30" 


2004-02-18 02:43 


2004-02-18 08:16 


17872 s 


30" 


2004-03-20 22:13 


2004-03-21 04:47 


20213 s 


22.5" 


2004-08-20 04:30 


2004-08-20 10:04 


18936 s 


30" 



that the source is a High Mass X-ray Binary (HMXB) at a dis- 
tance between 1 and 6 kpc, the mass donor being an early- 
type star, probably a sgB[e] star, surrounded by a dense and 
absorbing circumstellar material. The sgB[e] stars are massive, 
evolved, high-luminosity stars undergoing mass loss in a two- 
component wind (see, e.g., Hynes et al. 2002). 

The prominent, large equivalent width (EW) of the Fe 
K a , Fe K^, and Ni K a fluorescent emission lines, observed 
in the XMM-Newton spectrum of IGR J16318-4848, makes 
this object one of the most extreme examples of the new class 
of highly absorbed X-ray binaries discovered by INTEGRAL 
[see, e.g, Kuulkers, 2005 for a review]. 

A XMM-Newton observation campaign was carried out 
during 2004. At the same time, INTEGRAL was also observ- 
ing the source, collecting as well, data at energies larger than 10 
keV. Due to its prominent, large EW Fe and Ni emission lines, 
IGR J16318-4848 is an almost unique laboratory to study the 
physical condition and geometrical distribution of matter ac- 
creting on the compact object in a binary system. Thanks to 
this monitoring campaign we are able to compare the results of 
the radiation transfer simulations by Matt et al. (1999) and Matt 
(2002) to the spectra obtained through the coordinated obser- 
vations, in order to determine the time evolution of the follow- 
ing quantities characterizing the accretion flow: column density 
along the line of sight, covering fraction of the line-emitting 
matter, some constraints on the geometry, and ionization state 
of the matter. 

This paper is organized as follows: in section [2]information 
about the observations and the data reduction are given, while 
in sections [3]and |4]the results are presented and discussed. 

2. Observations and data reduction 

2.1. XMM-Newton 

The monitoring campaign consists of three observations spaced 
by periods of one month and five months, respectively. In this 
paper, we present a detailed study of these three observations 
along with the reanalysis of the XMM-Newton observation per- 
formed in February 2003. Table Q] summarizes the information 
on all XMM-Newton observation covered by this paper, to- 
gether with the EPIC-pn total exposure time and the extraction 
region used for each observation. The size of the extraction re- 
gion for each observation was selected such as to optimise the 
signal to noise ratio. 



In this paper, we only discuss the EPIC-pn data (Striider et 
al. 2001) as this instrument has the largest effective area above 

5 keV. We have analysed the EPIC-mos data and found that 
it does not contribute significantly to the results of this paper. 
Also, OM data was not considered since it is outside the scope 
of this paper. 

All data sets were analyzed using SAS v.6.1.0 (Gabriel et 
al. 2004). Following the standard procedure for XMM-Newton 
data reduction, we used the most recent calibration files avail- 
able in September 2004. 

The spectra were accumulated using single- and double- 
events (events where the charge in the CCD is collected in ei- 
ther one or two adjacent pixels). Pile-up effects were checked 
using the epaplot task and were found to be negligible in all 
observations. The background contribution was always small, 
except in the observation of March 2004, which was affected 
by high background radiation. Time intervals for which back- 
ground flaring does not improve the signal to noise in the EPIC- 
pn observations were excluded from the analysis following the 
method described in Piconcelli et al. (2004). TableQ]shows the 
exposures times after correction for background flaring activ- 
ity. 

All the spectra have been rebinned, in order to ensure that 
each background-subtracted spectral channel has at least 25 
counts, and that the instrumental energy resolution is not over- 
sampled by a factor larger than 3. 

2.2. INTEGRAL 

In Table [2] we present the observing log of the INTEGRAL 
(Winkler et al. 2003) observations which were performed si- 
multaneously with XMM-Newton. In the first two observations 
IGR J16318-4848 was the main target and the observations 
were done using a so-called 5x5 rectangular dither pattern. In 
the third observation the main target was the nearby (~ 56' 
away) source IGR J16320-4751 (Rodriguez et al. 2006) and 
the observation was done using a so-called hexagonal dither 
pattern. The 5x5 rectangular dither pattern consists of a square 
pattern around the nominal target location (1 source on-axis 
pointing, 24 off-source pointings, with a 2 degrees step), while 
the hexagonal dither pattern consists of a hexagonal pattern 
around the nominal target location (1 source on-axis pointing, 

6 off-source pointings, each 2 degrees apart). IGR J 163 18-4848 
was always in the fully coded field of view of both the IBIS and 
SPI instruments, where the instrumental response is optimal. 



Table 2. INTEGRAL monitoring campaign observations on 
IGR J16318-4848. The table shows the start time and end time 
of IBIS/ISGRI INTEGRAL observations on this source. 



Start Time - End Time (UT) 

2004-02-18 03:41 2004-02-18 17:29 
2004-03-20 20:15 2004-03-21 12:35 
2004-08-19 13:46 2004-08-20 22:41 
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Fig. 1. XMM-Newton EPIC-pn light curves in the 8-15 keV (upper panel) and 6-6.6 keV (middle panel) energy ranges. In the 
lower panel, the ratio of the two light curves is shown (bin time=1024 sees). 



In this paper we only consider data from IBIS/ISGRI 
(Lebrun et al. 2003 and Ubertini et al. 2003) and use only 
those data that overlap with the XMM-Newton observations. 
We do not consider the data from the IBIS/PICsIT, SPI, Jem-X 
or OMC instruments, since either the angular resolution is too 
high (SPI: 2.5°) and therefore the various sources in the Norma 
region of the Galaxy close to each other make a proper analysis 
not possible (e.g., Rodriguez et al. 2006), or IGR J16318-4848 
was too weak to be detected (IBIS/PICsIT, Jem-X, OMC). 

The INTEGRAL data reduction and the extraction of the 
events were performed using one of the latest versions of the 
Off-line Scientific Analysis (OSA; Courvoisier et al. 2003), i.e., 
v5.0. The data from the IBIS/ISGRI instrument were used to 
produce images in the 20^-0 and 40-80 keV energy ranges, 
with the aim of effectively identifying the most active sources 
in the field. From the results of this step, we produced a cata- 
logue of active sources which was given as an input for a sec- 
ond run producing images in the 20-60 and 60-200 ke V energy 
ranges; we increased the first range up to 60 keV to increase the 
signal to noise ratio in the lower energy band. In the latter we 
forced the extraction/cleaning of each of the catalogue sources 



in order to obtain the most reliable results for IGR J 163 18- 
4848 (see Goldwurm et al. 2003 for a detailed description of 
the IBIS analysis software). These results were used in the pro- 
duction of the energy spectra for IGR J16318-4848. Note that 
in all these processes, in order to be completely simultaneous 
with the XMM-Newton, we restricted the data reduction to the 
times that were strictly simultaneous with the XMM-Newton 
observations (see Table [T}. To do so we created good time in- 
terval fits files that were provided as an input to OSA; we also 
restricted the spectral energy range to 20-60 keV, since no sig- 
nificant emission was detected above these energies. 

3. Results 

3.1. Timing analysis 

In Fig. [Tj we show the IGR J16318-4848 light curves in the 
energy ranges, 8-15 keV (upper panel) and 6-6.6 keV (middle 
panel), together with the ratio of the two light curves (lower 
panel). We used a bin time of 1024 sees, although we also stud- 
ied several other bin times, such as, 256 and 512 sees, but they 
do not provide further improvement on the data quality. In this 
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Fig. 2. XMM-Newton count rate divided by the average count rate in each band. The rapid flux increments in the February 2003, 
February 2004 and August 2004 observations take places simultaneously in both bands. 



plot we can clearly see the February 2003 high energy flare that 
was studied by Matt & Guainazzi (2003) and variations by a 
factor 2 ~ 3 in < 5 x 10 3 s are visible also in the February 2004 
light curve, and a short flare in the August 2004 light curve. 

Only the February 2003 observation shows a significant 
variation in the ratio of the two light curves. This effect could 
indicate variations of the properties of the cold matter on times 
scales of about 10 4 s (Matt & Guainazzi, 2003). This variation 
in the continuum flux with respect to the flux in the Fe K a line 
is not present in the rest of the observations. 

In Fig. [2] we show the count rate in the 8-15 keV energy 
range, which is dominated by the continuum, and the count rate 
in the 6-6.6 keV energy range, which is dominated by the Fe 
K a line, divided by the average count rate in each band. During 
the February 2003, February 2004 and August 2004 observa- 
tions, a flux increment can be clearly seen. In all cases the Fe 
K a line follows almost immediately the continuum. This sug- 
gests that the emission lines are produced by the illumination of 
optically-thick matter on a small scale around the high-energy 
continuum source. These results are consistent with Walter et 
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band as a function of the 8-15 keV count rate. 
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Fig. 4. Spectra of the four IGR J 1 63 1 8-4848 XMM-Newton ob- 
servations. Despite the large flux variations between different 
observations, the spectral shapes follow the same trend. 

al. (2003). However, the variation of the continuum exhibits a 
larger dynamical range than the variation of the Fe K„ iron line 
(see Fig. [3}. This suggests that at least part of the line emission 
comes from regions, which do not respond on small timescales 
to the continuum. 

3.2. Spectral analysis 

In Fig. |4]all four IGR J16318-4848 spectra are shown. Despite 
the large flux variability, the spectral shapes are very similar 
one another. All the spectral analysis has been made using 
events above 5 keV because the source is highly absorbed be- 
low this energy. Detections at a level larger than 3 cr have been 
found only in the February 2004 and August 2004 observations 
(in Table|3]we shown the count rate (CR) in the 1-5 keV energy 
band). 
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Fig. 5. T - Nh ISO-^f 2 contour plots corresponding to IGR 
J16318-4848 XMM-Newton observations. 



Following Matt & Guainazzi (2003), the spectra of all avail- 
able XMM-Newton observations of IGR J 163 18-4848 have 
been fitted (using XSPEC v. 11. 3.1), alongside the simultane- 
ous INTEGRAL spectra whenever available, with the follow- 
ing baseline model: 

4 2 

F(E) = e- NH ^ E) [AE- T e EIE ' + ^ G t {E) + ^ S „,(£)] 



where cr(E) is the sum of the photoelectric and the Thompson 
cross-sections, Gj(E) are Gaussian profiles corresponding to 
the Fe K„ doublet, Fe and Ni K a fluorescent emission lines, 
and S CSi {E) is the Fe K ff i and Fe K a 2 Compton Shoulder re- 
spectively. The energy first component of the Fe K ff doublet 
has been left free to vary in the fit, while the energy centroid of 
the second component of the Iron line doublet was fixed to the 
energy of the first component multiplied by 1 .002 (Kranse et al. 
1979), as for neutral iron. The relative intensity of the second 
component of the Iron doublet was fixed to twice the intensity 
of the first component, as dictated by atomic physics. 

The last component of the baseline model was modeled 
using two box functions (as a first order approximation) cor- 
responding to the Compton Shoulder (CS) of Fe K ff doublet 
(e.g. Sunyaev & Churazov 1996, Matt 2002; Watanabe et al. 
2003). The first box function corresponding to the lowest en- 
ergy component of the iron doublet was extended in the closed 
interval [£i-0.16 keYZsi-O.Ol keV] as appropriate for the one- 
scattering CS (Matt 2002), where E\ is the best-fit centroid en- 
ergy for the Fe K ff i line, and 0.16 keV is the maximum recoil 
energy. We neglected the higher order shoulders, i.e. the ones 
dominated by two or more scatterings, because they are much 
fainter even for Compton-thick media, both in the reflection 
(Matt et al. 1991; George & Fabian 1991) and in the transmis- 
sion (Leahy & Creighton 1993) cases. Similarly, the box model 
for the second component was extended in the closed interval 
[£ 2 -0.16keV,£ 2 -0.01 keV]. 

In principle, one should include in the fit the Compton 
Shoulder for the other fluorescent lines as well. However, these 
lines are much fainter, and their Compton Shoulder are basi- 
cally undetectable in the EPIC-pn spectra. Therefore, for the 
sake of simplicity they will not be considered in what follows. 
For the same reason, the Fe and Ni K a lines were modeled 
as single Gaussians. 

Regarding the absorbing matter, the abundance (normal- 
ized to the solar value) of elements with Z > 26, (Zp e here- 
after) was allowed to vary independently of the abundance 
of the lighter elements. The relative iron abundance values 
(Anders E. & Grevesse N.) obtained for each observation were: 
0.81 + 0.03, 0.77 ±g;[g, 0.89+°^ and 0.86±°-{* for the February 
2003, February 2004, March 2004 and August 2004 observa- 
tions, respectively. These values are compatible within the er- 
rors and therefore we calculated the weighted mean, which re- 
sulted to be 0.81 + 0.02, and used it in all following spectral 
fits. 

XMM-Newton and INTEGRAL spectra were fitted si- 
multaneously applying the relative normalization between the 
EPIC-pn and ISGRI data given by Kirsch et al. (2005). 
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Table 3. Best-fit parameters and results of the IGR J16318-4848 monitoring campaign using the baseline model. 



Parameter 



February 2003 February 2004 March 2004 August 2004 



1 c 4 +U.02 

x ""-0.4 

35 ± 1 



2.12 ±0.03 
6.39 ±0.01 
1.71 ±0.02 
9.5 ± 0.3 
7.05 ± 0.01 
0.20 ± 0.02 

9 9+0.3 
^-0.2 

7 45+0.O6 
' -o.oi 
05 +0 02 

7 +0 1 
"•'-o.i 

0.05 ± 0.03 
106.9/92 

120.9/93 

9 oq+o.03 

• £ " J -0.06 

1.3 ±0.6 



'~ 4 cm- 2 s M 



r 

E c [keV] 

A [lO" 2 cm" 2 s -1 ] 
HR[|5|] [keV] 
N H [10 24 cm" 2 ] 
E x (Fe K ffl ) [keV] 
/i(FeK„)[10 ' 
EW (Fe K ff ) [eV] 
E 2 (Fe K» [keV] 
h" 

EW (Fe K/j) [eV] 
£ 3 (NiK tt ) [keV] 

EW (Ni K a ) [eV] 
j a 

2 i 

X 2 /v (Without CS) 

F [10-" erg cm" 2 s -1 ] 

CR (1 - 5keV) [10- 3 cts s" 1 ] 



1.46 ±0.03 
< 16 
27.0!™ 
0.40 ±0.14 

1 Q9+0.05 

6 43+ - 01 

u - H -'-0.02 



0.93 : 
8.0: 



0.04 
0.4 



7 10+ 003 

'• lu -0.02 

24+ 03 

9 a+0.3 

^-0.5 

7 45 +ai 
' -o.oi 

05+ 03 
5 +04 

-0 2 

0.25 ± 0".06 
98.9/97 
153.0/98 
2.11 ±0.08 
12.9 ± 1.4 



1.35 ±0.04 

< 180 
9.95 ±0.01 
0.46 ± 0.04 

2 16 +007 
6.44 ± 0.01 



0.8!»J 
<2.5 



1.41 ±0.03 
< 30 
4 8 +01 

^•°-0.2 

0.44 ±0.10 

1 QQ+0.04 
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6 42 +a05 

-0.04 





65 +0 06 

"• UJ -0.03 
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7 c+0.4 
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7 c+0.4 
' -0.7 


03 +0 02 


07+ 03 
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5 +04 
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9 9+0.8 
^-1.3 


24 +0 13 

"•^-0.06 


19 +008 


64.9/51 


111.7/83 



114.85/52 135.85/84 



1 2 +0 - 3 
2.10 ±0.40 



"ratio intensity against the Fe K a 



We studied the width of the emission lines in order to avoid 
any possible numerical problem related to unresolved lines. 
Firstly, we fit all the spectra with the line width fixed to 10~ 4 
eV. We then allowed the line width of the emission lines to vary, 
but no improvements in the x 1 was found. With these condi- 
tions, we obtained very tight upper limits for the width of the 
emission lines, smaller than 10 eV. The only exception is the 
August 2004 observation, when the width is basically uncon- 
strained. We then fixed the line width to 10~ 4 eV for all the 
observations. 

Despite the strong coupling between V and Nh that was 
found in the February 2003 observations (Matt & Guainazzi, 
2003), we have been able to find a stable spectral index for 
all the observations, with the exception of August 2004 (see 
Fig. |5). This has been posible due to the large energy cover- 
age introduced by the simultaneous INTEGRAL data; the bet- 
ter constraint on the metal abundance and the inclusion in our 
model of the Fe K a doublet (see Pig. [SJ. 

In Table [3] we summarize the best-fit parameters of the 
baseline model and in Fig.|6]we present the spectra and residu- 
als of the data against the baseline model. All the errors shown 
in this paper correspond to a 90% confidence level for one in- 
teresting parameter. The values of the equivalent widths in the 
table were computed with respect to the unabsorbed contin- 
uum. 

In all 2004 observations, the relative intensity of the 
Compton Shoulder is larger than in the 2003 observation. 
Interestingly, for these observations the centroid of the Fe 
Kq,i line is significantly bluer than expected for neutral iron 
(6.40 keV). This seems to suggest significant ionization of iron 
(Fe > XIV, Kallman et al. 2004). However, the strong iron 
Fe K^g line implies that iron is less ionized than Fe xi or so 
(Molendi et al, 2003). The apparent blueshift may then be an 
instrumental effect, or an artifact due to either the oversimpli- 



fied CS modeling or to the assumption of a narrow line core. 
We recall that the systematic uncertainties on the gain of the 
EPIC-pn camera at 6 keV is 10 eV (Kirsch 2006). 

In Fig. [7] we show the residuals of the Compton-Shoulder 
feature for all the observations. We have also included in 
the Table [3] the x 2 values of our baseline model without the 
Compton Shoulder. In all the observations the improvement 
adding the CS is statistically significant. 

A comparison between the ratio of all INTEGRAL spec- 
tra and their mean spectrum was performed, showing no trend. 
This stable spectral shape agrees with the accreting pulsar ob- 
jects. 

Interestingly, Matt & Guainazzi (2003) report soft excess 
emission above the extrapolation of the baseline model at ener- 
gies E < 5 keV. We have reanalyzed the images obtained dur- 
ing the monitoring campaign and discovered a serendipitous 
source at ^30" from the IGR J16318-4848 centroid position. 
This serendipitous source is totally responsible for the reported 
soft excess. In the other three observations, this soft source falls 
into the CCD gaps. Still, we detect soft X-ray emission at a sig- 
nificance level larger than 3 <x in two out of the four observa- 
tions. No spectral information can be derived on this compo- 
nent from the available data. 



4. Discussion 

We have presented results from a coordinated XMM- 
/Vewfon/INTEGRAL monitoring campaign on IGR J 163 18- 
4848. The source was always clearly detected, albeit with large 
amplitude (factor of 3) flux variations, suggesting (along with 
its rediscovery in the ASCA archive) that it is a persistent 
source which escaped attention before 2003 only due to the 
strong absorption. Persistence of the X-ray emission is a char- 
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Fig. 6. Spectra (upper panels) and residuals in units of standard deviations (lower panels) when the baseline model is applied to 
the IGR J16138-4848 observations discussed in this paper. 
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Fig. 7. Residuals of the Compton-Shoulder feature for all of 
the IGR J16318-4848 observations. 



acteristic of Be systems with wide orbits ( hundreds of days, 
Negueruela 1998). 

The properties of the circumstellar matter clearly changed 
from one observation to the other, the column density (as de- 
rived from absorption) for instance being significantly lower 



(but still Compton-thick) in the last observation. As a further 
probe of insight on the circumstellar matter, we used as a diag- 
nostic the Compton Shoulder, which depends on the covering 
fraction of the matter and on its average column density (Matt 
2002, Matt & Guainazzi 2003), while the absorption of course 
gives information on the line-of-sight (l.o.s.) column only. 

In Fig. [8] we show the results of the Fe K a Compton 
Shoulder normalized intensity, fcs (the ratio between the CS 
and the line core intensities) against the EW of the Fe K„. 
Comparing our results with the predictions of radiative trans- 
fer simulations done in Matt (2002), a covering fraction of the 
obscuring matter in the range 0. 1-0.3, and average column den- 
sities in the range 2 x 10 22 to 5 x 10 23 cirT 2 are derived. 

It is worth noting that the average column density is always 
much smaller (1 to 2 orders of magnitude) than the l.o.s. one. 
The fact that the two values are largely different implies that 
the circumstellar matter is highly inhomogeneous. As the av- 
erage one is systematically smaller that of the l.o.s. it strongly 
suggests that we are observing the source from a privileged di- 
rection. This, together with the relatively small covering factor 
(which rules out the hypothesis that the compact object (a neu- 
tron star according to Filliatre & Chaty 2004, c.f., Negeruela 
1998) is embedded in the wind from the giant companion (a 
sgB[e] star according to Filliatre & Chaty 2004), not surpris- 
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Fig. 8. Fe K [V fluorescent line EW against the Compton 
Shoulder normalized intensity (fcs)- The lines represent the 
predictions of radiative transfer simulations after Matt (2002) 
for different values of the gas covering fraction, C/. The filled 
squares represent the four XMM-Newton measurements dis- 
cussed in this paper. The empty squares indicate the position 
along the lines corresponding to the following values of col- 
umn density of the obscuring gas: N H = 10 22 , 2x 10 22 , 5 x 10 22 , 
10 23 , 2 x 10 23 , and 5 x 10 23 cirT 2 . 
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ingly if the system has indeed a wide orbit, as suggested by its 
persistency), can be interpreted in terms of a flattened distri- 
bution of matter (quite plausible in a binary system, especially 
if with wide orbit), with the column density decreasing from 
the equatorial plane, and the line of sight lying close to that 
plane. Even if this scenario is by no means unique, it seems to 
be the simplest and most natural. The variations in the cover- 
ing factor and average column density may then be the witness 
of variations in the mass transfer rate, possibly due to a highly 
eccentric orbit of the system. 
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